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We study the behaviour and consequences of cosmic string networks in contracting universes. They 
approximately behave during the collapse phase as a radiation fluids. Scaling solutions describing 
this are derived and tested against high-resolution numerical simulations. A string network in a 
contracting universe, together with the gravitational radiation it generates, can affect the dynamics 
of the universe both locally and globally, and be an important source of radiation, entropy and 
inhomogeneity. We discuss possible implications for bouncing and cyclic models. 


INTRODUCTION 

Cosmological scenarios involving oscillating or cyclic 
universes have been know for a long time [jj] . Recent in¬ 
terest has been associated with a cyclic extension of the 
ekpyrotic scenario (Q. A related result was the realiza¬ 
tion &!§ that the presence of a scalar field seems to be 
necessary to make cosmological scenarios with a bounce 
observationally realistic. And if scalar fields are present, 
we should contemplate the possibility of topological de¬ 
fect formation [Q. Here we study cosmic string evolution 
in a collapsing universe, and discuss some implications of 
the presence of cosmic strings for bouncing universes. In 
a bouncing universe scenario the properties of the uni¬ 
verse in the expanding phase depend on physics happen¬ 
ing in a previous collapsing phase. Hence, if defects do 
exist in these models, it is crucial to understand their 
evolution and consequences in both the expanding and 
collapsing phases. In particular, we expect that cosmic 
strings will become ultra-relativistic, behaving approxi¬ 
mately like a radiation fluid. This means that a cosmic 
string network, both directly and through the gravita¬ 
tional radiation emitted by its loops, will soon become a 
significant source of entropy (and inhomogeneity), mak¬ 
ing it a further problem for cyclic universes if a suitable 
and efficient mechanism for diluting the entropy is not 
available. A more detailed analysis can be found in (7j|. 


COSMIC STRING EVOLUTION 

The world history of a cosmic string can be represented 
by a two-dimensional surface in space-time, obeying the 
usual Goto-Nambu action, from which it is easy Q to 
derive the microscopic string equations of motion. Con¬ 
sider for a start the evolution of a circular cosmic string 
loop in a cyclic universe. Simple analytic arguments show 
that a loop whose initial radius is much smaller than the 


Hubble radius will oscillate freely with a constant invari¬ 
ant loop radius and an average velocity v = I/a/ 2- (Note 
that we are assuming units in which c = K = 1.) On 
the other hand, once the collapse phase begins, we will 
eventually get to a stage in which the physical loop ra¬ 
dius becomes comparable to the Hubble radius ar ~ H~ 1 
and then gets above it. In this regime the loop velocity 
is typically driven towards unity v —> 1 and it straight¬ 
forward to show that the invariant loop length (which is 
proportional to the energy of the loop) grows as R oc a -1 
and the Lorentz factor as 7 oc a~ 2 . Despite its growing 
energy R , the actual physical loop radius ar = R/j —> a, 
so the loop shrinks with the scale factor and inexorably 
follows the collapse into final big crunch singularity. 


Importantly, this relativistic final state for a loop in a 
collapsing universe is generic and quite different to the 
initial condition usually assumed for super-horizon loops 
in the expanding phase. There, loops begin with a van¬ 
ishing velocity which only becomes significant when the 
loop falls below the Hubble radius. Such evolution can¬ 
not be reproduced in reverse during the collapsing phase 
without fine-tuning the velocity as the loop crosses out¬ 
side the Hubble radius. This simple fact introduces a 
fundamental time asymmetry for string evolution in a 
cyclic universe (and for all other defects). The analytic 
expectations for our circular loop solution can easily be 
confirmed by a numerical study |7|. 


Two complementary approaches are available to study 
the evolution of a cosmic string network: one can re¬ 
sort to large numerical simulations J. I 0 , or one 
can develop analytic tools [jfl], 0 0 which provide an 
averaged description of the basic properties of the net¬ 
work. We shall use the best motivated of these analytic 
models, the velocity-dependent one-scale (VOS) model 
@00 [Tc| . The VOS model describes the string dy¬ 
namics in terms of two ‘thermodynamical’ parameters: 
the string RMS velocity, Vqo, and a single length scale— 
the string network is thus assumed to be a Brownian 
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random walk on large enough scales, with a correlation 
length L. Hence one can simply relate it with the energy 
density in long strings as p x = p/L 2 , where p is the 
string mass per unit length. Note that the commonly 
used ‘correlation length’ L is really a measure of the in¬ 
variant string length or energy, rather than the typical 
curvature radius of the strings. By including the appro¬ 
priate Lorentz factor 7^ = (1 — v j^,) -1 / 2 f° r the long 
strings, we can denote the physical correlation length by 
iphys = £700 • With this assumption the VOS model has 
one phenomenological parameter c, commonly called the 
loop chopping efficiency, which describes the rate of en¬ 
ergy transfer from the long-string network to loops. The 
evolution equations then take the following form 0 ® 

2—= 2(1 + v^HL + cVoo + STGpv^, (1) 

dVoO 2 / k{ V Oo) r> jj /V,\ 

-*- = < 1 ' (2 > 

The final term in the evolution equation for the corre¬ 
lation length describes the effect of gravitational back- 
reaction. We are not including in either equation addi¬ 
tional terms arising from friction due to particle scatter¬ 
ing p6| , which could conceivably be important during 
the final stages of collapse. We shall return to this point 
below. Here fc(i’oo) is the momentum parameter, which is 
thoroughly discussed in [jn|. Notice that this is positive 
for 0 < v 2 < 1/2 and negative for 1/2 < v 2 < 1. 

When the contraction phase starts and the Hubble pa¬ 
rameter becomes negative the velocity will tend to in¬ 
crease: as in the simple case of the circular loop, the 
string velocity will gradually tend towards unity. In this 
approximation, and neglecting for the moment the loop 
production and gravitational back-reaction terms, the 
evolution equation for the correlation length easily yields 
L oc a 2 . Note that this is the same overall scaling law for 
the string network as that in a radiation-dominated ex¬ 
panding universe; the string network effectively behaves 
like a radiation component. In terms of the physical cor¬ 
relation length, Lphys oc a, as if strings were being con¬ 
formally contracted (except for their rapidly growing ve¬ 
locities). However, there are several factors that must be 
considered which complicate this simple state of affairs. 
First, there is the issue of loop production. Under the 
above assumptions on velocity, but putting the loop pro¬ 
duction term back in the correlation length equation, one 
finds the following approximate solution in the radiation 
and matter eras, 

d-'rad = ( Lmax — In CL 1 d (3) 
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where L max is the string correlation length at the time 
of maximal size of the universe, and the scale factor at 


that time was chosen to be unity (so the logarithm term 
in the first case is positive). Hence if c remains constant 
(or is slowly varying), asymptotically the scale factor de¬ 
pendent terms will dominate, so that L oc a 2 In a in the 
radiation era, and L oc a 3 / 2 in the matter era. The latter 
is also the scaling law for the correlation length in the 
matter-dominated, expanding universe. This highlights 
the different roles played by loop production in the scal¬ 
ing behaviour of a cosmic string network in the radiation 
and matter eras 00- A strong argument can be made, 
however, for a relativistic correction to the loop produc¬ 
tion term. In the simplest form of the VOS model there 
is an identification between the correlation length, L , and 
the physical distance T p hy S which a string segment is ex¬ 
pected to travel before encountering another segment of 
the same size forming a loop in the process. However, 
taking into account the Lorentz factor in the physical 
correlation length, one would expect 0 that c oc 7oo 1//2 
thus driving c rapidly towards zero and asymptotically 
yielding our simple solution L oc a 2 both in the radia¬ 
tion and matter eras. Of course, during re-collapse we 
expect that c will depend on a number of other proper¬ 
ties of the string network such as the enhanced build-up 
of small scale structure due to the contraction. Even¬ 
tually, however, the Hubble radius will fall below even 
the length scale of wiggles on the string after which our 
asymptotic solution L oc a 2 should be valid. In what 
follows, we shall consider the two well-motivated cases, 
first, c = const. ^ 0 initially and, secondly, c = 0 the 
probable asymptotic case. Further supporting evidence 
for this behaviour is discussed in [f 7 l|. 

Returning to our analytic solutions for the constant 
loop production case 1) and (||), we can use the velocity 
equation to find an approximate, implicit solution 


i 2 4 

1 — v oc a exp 


2 k(v) f a l t x dci 


( 5 ) 


A J L(a) J ’ 

where A = 1 in the radiation era and A = 2 in the matter 
era. Substituting (j|-[fj) one respectively obtains 

\4 k{v)/c 


1 - Vrad a (-Ina) 

a 4 +W. 


( 6 ) 

( 7 ) 


Hence in the limit where v —> 1 and therefore k —> 0 
the asymptotic solution would have the form 7 -2 oc 
(1 — v 2 ) oc a 4 . The momentum corrections, which phe¬ 
nomenologically account for the existence of small scale 
structures on the strings, imply that convergence will be 
slower than this. These solutions will still hold when one 
includes the gravitational back-reaction term jlGj . As a 
final caveat, it is also worth emphasizing that the VOS 
model assumes that the long string network has a Brown¬ 
ian distribution on large enough scales, which may not be 
a realistic approximation in a closed, collapsing universe. 
This point clearly deserves further investigation. 
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As a test to the above solutions, we have performed 
a number of very high resolution Goto-Nambu simula¬ 
tions on the COSMOS supercomputer, using a modified 
version of the Allen-Shellard string code |J. Further nu¬ 
merical details can be found in [|7J. Our results are con¬ 
sistent with the existence of an attractor solution of the 
type described above. The result of two such simulations, 
for universes filled with radiation and matter, is shown 
in Fig. |l|. During the expanding phase we confirm the 


usual linear scaling regimes in 

the radiation and matter 

eras, respectively 



Gexp,rad OC t CX O , 

Vqq = const. 

(8) 

T + 3/2 

■^exp^mat OC Z OC Cl , 

Voo = const. 

(9) 
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Once the contraction starts, these are modified: the 
velocity starts increasing, and the scaling of the cor¬ 
relation length with the scale factor also drops, being 
approximately constant to begin with, and then rising 
slowly. One can identify a transient scaling phase, valid 
in the period p ~ 1.0 — 1.4, where one approximately 
has Ltrans oc a in the radiation-dominated case, and 
Ltrans oc a 5 / 4 in the matter era. Unfortunately, the ex¬ 
tremely demanding requirements in terms of resolution of 
the simulation do not currently allow us to run simula¬ 
tions with longer dynamic range to establish beyond rea¬ 
sonable doubt whether this scaling law approaches j3 = 2, 
as predicted above. However, there are strong indications 
that the networks are evolving towards this asymptotic 
regime, as shown by the relatively rapid climb of the ex¬ 
ponent in Fig. jjj It is clearly noticeable that the velocity 
rises much faster in the matter era than in the radia¬ 
tion era. It is also interesting to point out that during 
the collapse phase the loop and long string velocities are 
noticeably different, and this difference (which is more 
significant in the radiation than in the matter case) in¬ 
creases with time. The plot also shows an apparent dif¬ 
ference in the expanding phase, but this is not significant: 
the initial lattice conditions tend to give different veloc¬ 
ities to small loops than to long strings and they start 
evolving, and this difference is gradually erased. 

We also notice that the network keeps chopping off 
loops throughout the simulation, and that there is a dra¬ 
matic increase in the small scale structure of the net¬ 
work, particularly at later times. Visually, the string 
network develops large numbers of ‘knots’, highly convo¬ 
luted strings regions where the wiggly long strings have 
collapsed inhomogeneously. These small scale features 
have proved to be difficult to evolve numerically, and 
this in fact turns out to be the main limiting factor at 
present preventing us from running the simulations closer 
to the big crunch. A comparison of our solutions (|3|- 
and (](]-[7|) with the numerical simulations described 
above, produces a very good matching—see Jt[ . Finally 
it is also worth keeping in mind that any discussion of 


FIG. 1: Cosmic string network properties in the expansion 
and contraction phases, as a function of conformal time (such 
that maximal expansion occurs at rj = 1 and r] crurlc h = 2). 
Solid is a simulation in the matter era, dashed is radiation. 
Plotted are the re-scaled correlation length L/t (top left), the 
scaling law of the correlation length L oc a 13 (top right), the 
velocity of the long string network (middle left), the ratio of 
the loop and long string velocities (middle right), the fraction 
of the string energy in the form of long strings (bottom left) 
and the total number of loops (bottom right). 

the evolution of a cosmic string network with the present 
formalism is only applicable while one is well below the 
Hagedorn temperature, at which the strings would ‘dis¬ 
solve’ in a reverse phase transition. Discussions of asymp¬ 
totic regimes should be taken with some caution, since 
a cosmic string network will only survive the bounce in¬ 
tact if this happens before the Hagedorn temperature is 
reached. 


DISCUSSION: COSMOLOGICAL 
CONSEQUENCES 

The overall density of strings remains constant rela¬ 
tive to the background density p in both radiation and 
matter eras, Poo/p = crGp, with er r ss 400 and er m « 60 
respectively || [|. During curvature domination or ac¬ 
celerated expansion, the string density grows relative to 
the other matter as Poo/Pm oc a. For GUT-scale strings 
with Gp ~ 10~ 6 this gives the interesting conclusion that 
today strings have a comparable energy density to the 
CMBR. However, a realistic cyclic model will continue 
to expand well beyond foi so the string density at max¬ 
imum expansion will end up being much greater than 
the radiation density. In addition, the gravitational (or 
other) radiation produced through the continuous decay 
of the string network evolves as p gr oc a~ 4 . It might ap¬ 
pear that this contribution would become negligible dur¬ 
ing the matter era but in each Hubble time the strings 
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lose about half their energy into gravitational radiation, 
so this background always remains comparable to the 
string density p gr ~ p^. Now consider the collapsing 
phase in which the string network, like the gravitational 
waves they have produced, begins to behave like radia¬ 
tion. Globally, the density of both the strings and the 
gravitational waves will grow as a ~ 4 and, together with 
any other radiation components, they will eventually 
dominate over any nonrelativistic matter. In a realistic 
cyclic model, sufficiently massive strings and their decay 
products will have a greater density than the microwave 
and neutrino backgrounds. As the universe contracts, 
it will eventually reach a state in which the relativistic 
string network and/or their gravitational waves dominate 
the global dynamics of the universe! This would lead to a 
dramatically different universe after it emerges from the 
next bounce. Even lighter strings, which do not dominate 
the universe, would end up with a much greater density 
in the collapsing phase than they had previously during 
expansion. If the universe went through a bounce, the 
energy density in the cosmic strings and gravitational ra¬ 
diation produced by the network would be much greater 
after the bounce than before it. Bounds on the string 
mass per unit length may be severely modified, in ad¬ 
dition to more general constraints on extra relativistic 
fluids JTt| . 

Furthermore, unlike the uniform CMB, the energy den¬ 
sity in both cosmic strings and gravitational radiation 
will be very inhomogeneous. In the collapsing regime, an 
increasingly small fraction of Hubble regions will have a 
string passing through them. Those that do will become 
string dominated since the string energy density in those 
regions will approximately evolve as Poo/p cx 700 oc a~ 2 , 
up to the corrections described above. For these regions 
the assumption of a FRW background will cease to be 
valid at late times, and the defects can make the uni¬ 
verse very inhomogeneous fis)] oranisotropic ^J. Even 
Hubble regions without strings will have large fluctua¬ 
tions in their gravitational radiation content. For suffi¬ 
ciently massive strings, both of these effects can survive 
the bounce to create large inhomogeneities in the next 
cycle. 

A possible caveat to these solutions is dynamical fric¬ 
tion (which we have neglected so far). In the 7 >> 1 
limit, one can estimate [Q that strings loose all their 
momentum due to this effects in one Hubble time when 
Gp ~ 7 -2 (radiation era result) or Gp ~ 7 _1 (mat¬ 
ter era result). However, this assumes a homogeneous 
and isotropic background, so need not apply in our case. 
Moreover, the fact that a significant amount of momen¬ 
tum will be transferred from the strings to the back¬ 
ground will in itself add to the anisotropies which natu¬ 
rally occur in our model. This will be developed further 
elewehere 0 . 

We conclude that a cosmic string network will be a sig¬ 
nificant source of radiation, entropy and inhonrogeneity 


which may be problematic in the cyclic context. Some of 
the results described in this paper should also be valid 
for other topological defects, in particular domain walls. 
Conversely, if direct evidence is found for the presence of 
topological defects in the early universe, their existence 
alone will impose constraints on the existence and char¬ 
acteristics of previous phases of cosmological collapse. 
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